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Abstract 

Background: Small airways are regarded as the elective anatomic site of obstruction in most chronic airway 
diseases. Expiratory computed tomography (CT) is increasingly used to assess obstruction at this level but there is 
no consensus regarding the best quantification method. We aimed to evaluate software-assisted CT quantification 
of air trapping for assessing small airway obstruction and determine which CT criteria better predict small airway 
obstruction on single breath nitrogen test (SBNT). 

Methods: Eighty-nine healthy volunteers age from 60 to 90 years old, underwent spirometrically-gated inspiratory 
(I) and expiratory (E) CT and pulmonary function tests (PFTs) using SBNT, performed on the same day. Air trapping 
was estimated using dedicated software measuring on inspiratory and expiratory CT low attenuation area (LAA) 
lung proportion and mean lung density (MLD). CT indexes were compared to SBNT results using the Spearman 
correlation coefficient and hierarchical dendrogram analysis. In addition, receiver operating characteristic (ROC) 
curve analysis was performed to determine the optimal CT air-trapping criterion. 

Results: 43 of 89 subjects (48,3%) had dN2 value above the threshold defining small airway obstruction (i.e. 
2.5% N2/I). Expiratory to inspiratory MLD ratio (r = 0.40) and LAA for the range -850 -1024 HU (r = 0.29) and for the 
range -850 -910 HU (r = 0.37) were positively correlated with SBNT results. E/l mld was the most suitable criterion 
for its expression. Expiratory to inspiratory MLD ratio (E/l mld) showed the highest AUC value (0.733) for small airway 
obstruction assessment. 

Conclusion: Among all CT criteria, all correlating with small airway obstruction on SBNT, E/l mld was the most 
suitable criterion for its expression in asymptomatic subjects with mild small airway obstruction 

Trial registration: Registered at Cli nicaltrials.gov, identifier: NCTOl 230879. 
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Background 

Chronic airway diseases are obstructive lung disorders oc- 
curring more and more frequently, therefore becoming a 
major public health burden worldwide [1]. Severity assess- 
ment and management of these diseases are defined ac- 
cording to clinical examination and routine pulmonary 
function tests (PFTs). However, these measurements are 
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not strong enough to accurately discriminate patients re- 
garding clinical outcomes [2]. Chronic obstructive pul- 
monary disease (COPD) and asthma are the most 
common causes of such chronic airflow limitation. From a 
pathophysiological point of view, small airways (defined 
from 2 mm in internal diameter and downward) are 
regarded as the underlying elective anatomic level of air- 
way obstruction for both disorders [3,4], but are inad- 
equately investigated using conventional PFTs, namely 
"the lungs quiet zone" [5]. Specific structural changes at 
this level can be assessed using tools such as the single- 
breath or multiple-breath nitrogen washout test (SBNT, 
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MBNT respectively) [6]. SENT and MBNT are time con- 
suming and the routine uses are not widely available 
owing to limited access to the equipment. Due to the fact 
that expiratory acquisition allows indirect evaluation of 
bronchiolar involvement [7], computed tomography (CT) 
has been increasingly put forward as an appropriate non- 
invasive tool for refinement in the classification and treat- 
ment monitoring of COPD and asthma [8-14]. The appli- 
cation of an increasingly wide range of technological tools 
also allows post-processing by segmentation software and 
thus quantifies air trapping objectively. Both a decrease in 
mean lung density and the percentage of low attenuation 
area on expiratory CT have been used in various study 
and correlated with disease severity in COPD and asthma 
[15-17]. Some authors have also suggested the role of 
paired expiratory to inspiratory ratio or difference as a 
small airway marker [18-20]. However, the choice of one 
of these criteria may affect the results and there is, there- 
fore, a crucial need for standardization in CT air trapping 
expression. 

This prospective study was designed firstly to explore 
the validity of air trapping software-assisted CT quantifi- 
cation for the assessment of small airway obstruction by 
using SBNT as standard of reference and, secondly, to 
assess which was the most accurate criterion. 

Methods 

Subjects and eligibility criteria 

Between August 2009 and April 2012, we prospectively 
conducted a cross-sectional study on lung aging. This 
specific population was chosen to select a range of sub- 
tle small airway obstruction. Eligible participants were 
asymptomatic, non-smokers for at least 20 years with a 
cumulative history <10 pack-years of tobacco use and 
without a past history of lung disease. Subjects were 
asked to be in good mental and physical health as 
assessed by medical interview and physical activity using 
the Voorips score questionnaire [21]. All the final study 
participants had normal range spirometry, specifically 
regarding forced expiratory volume. They were informed 
of the aims of the study and gave their informed consent 
for both pulmonary function testing and chest CT. 

This study received the approval of the local research 
ethics committee (CPP sud Mediterranee IV) and the 
agreement of the French Health Products Safety Agency 
(ANSM) before the start of the research. 

Scanning techniques 

All CT examinations were performed using a 64-de- 
tector row CT scanner (LightSpeed VCT; GE Healthcare, 
Waukesha, USA) without the administration of contrast 
material. Patients were placed in the supine position. 
The entire chest from apex to posterior recesses was in- 
cluded in the cranio-caudal direction according to the 



following protocol: tube voltage: 120 kV, automatic tube 
current modulation with maximal current limited to 300 
mAs, collimation: 64 x 0.6 mm, increment: 0.9, recon- 
structed slice thickness: 1.25 mm, tube rotation: 0.5 s, 
acquisition field of view ranged from 320 to 380 mm de- 
pending on the patients body habitus. Each chest CT 
examination was reconstructed using a standard filter. 

Images were acquired using spirometric gating (Win- 
spiroPRO; Medical International Research Waukesha, 
USA) to monitor lung volumes and determine respira- 
tory endpoints. Patients were trained to breathe into a 
mouthpiece connected to the handheld spirometer after 
calibration of the spirometer according to manufacturer 
specifications. Breath-hold was obtained at least 3 times 
in inspiration and expiration using the spirometric 
equipment just before CT acquisition in supine position. 
Image sets were then acquired during single-breath-hold 
helical scanning between 90 and 100% of previously re- 
corded slow vital capacity in supine position. Four add- 
itional slices were acquired in the same conditions at 
end expiration: one through the upper lobes, one 
through the tracheal carina and two between the tra- 
cheal carina and the diaphragm. All DICOM images 
were anonymously archived and transferred to a dedi- 
cated computer for post-treatment. 

Quantitative assessment of air trapping by CT 

De-identified data was processed independently and 
blindly by a senior radiologist with 8 years of experience 
interpreting chest CT images (S.B.) and also by a chest 
radiology fellow (F.K.). Lung parenchyma was automatic- 
ally isolated from the chest wall, mediastinum and air 
contained within the segmental bronchi (Additional file 
1: Figure SI) and then analyzed using threshold tech- 
niques and histogram computation in order to quantif)^ 
lung air trapping objectively and to provide reproducible 
data with the use of commercially available software 
(Myrian, Intrasense, Montpellier, France). Measurements 
were done for the four levels described above combining 
the values found for both sides at each level. Air trap- 
ping was defined in expiratory CT as lung regions that 
failed to increase in attenuation or decrease in volume 
in a normal fashion compared with findings on the ini- 
tial inspiratory images [22]. Given the absence of a sin- 
gle, validated method of CT air trapping quantification, 
the inspiratory and expiratory chest CT scans were 
scored in 3 ways, most frequently used in the literature 
[23]. 

The percentage of lung parenchyma that fell within a 
range of low attenuating area (LAA%) was considered 
representative of air trapping: voxels between -850 
and -1024 HU, called exp -850-1024 on expiratory 
slices were isolated. A second threshold was also applied 
between -850 and -910 HU (exp -850-910) to avoid 
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low-density values due to emphysematous or cystic 
lesions. Finally, mean lung density was automatically 
measured. 

Expiratory to inspiratory ratio and the difference be- 
tween these values were calculated. 4 expiratory images 
were matched visually by using anatomical landmarks 
with 4 inspiratory images with the same slice thickness 
from the inspiratory images. E/I, E-I and (E-I)/I were 
expressed using the two thresholds described above for 
LAA% and for MLD. 

Single breath nitrogen test (SBNT) 

SENT was performed the same day, just before CT scan- 
ning. This test was used to look specifically for small 
airway abnormalities. This test was performed as previ- 
ously reported by Bourdin and colleagues [24], using the 
Vmax apparatus (Vmax; Sensomedics, Yorba Linda, CA, 
USA) under the supervision of an expert in clinical 
physiology (M.H.), blinded to the CT results. Subjects 
were asked to breathe slowly and deeply to reach re- 
sidual volume. They were then asked to inhale as deeply 
as possible to obtain total lung volume capacity at which 
point the breathing valve started the delivery of 100% 
02 gas. Subjects were then asked to exhale immediately 
at a 0.3-0.5 1/s flow giving visual feedback on the com- 
puter screen to allow for the procedure to be followed 
correctly. The fractional expiratory nitrogen concentra- 
tion was plotted against the expired volume. This helped 
determine the closing volume and the slope of this rela- 
tionship when it reached the plateau known as the phase 
III nitrogen slope (delta N2 or dN2). See the Additional 
file 2: Figure S2 for the single breath nitrogen washout 
test curve. 

Statistical analysis 

Quantitative variables were expressed as means or 
medians and compared using the Student t test or the 
Wilcoxon test. Qualitative variables were expressed as 
numbers (percentages) and compared using the Chi 2 
test or the Fisher test as appropriate. 

Non-parametric Spearman test was carried out to as- 
sess the relationship between functional and CT param- 
eters. Distance and similarity between dN2 and CT air 
trapping criteria were then analyzed using dendrogram 
analysis based on the agglomerative hierarchical cluster 
tree. The best CT parameter was determined using re- 
ceiver operating characteristic (ROC) curve analysis 
using standard dN2 cut-off = 2.5% N2/1. To test the diag- 
nostic consistency of air trapping quantification for the 
two readers, inter-observer agreement was calculated 
using intra-class correlation coefficient [25]. A p-value 
of < 0.05 was considered statistically significant. The stat- 
istical analysis was performed using statistical software 



SAS 9.3 (SAS Institute; Cary, NC), R 2.14.1 and SPSS 
17.0 for Windows (SPSS, Chicago, 111). 

Results 

Among the 101 eligible subjects, a total of 89 performed 
adequate spirometric maneuvers for both PFTs and CT 
acquisitions. Software-assisted inspiratory and expira- 
tory lung segmentation was achievable for all these sub- 
jects. Description of the study population is provided in 
Table 1. 



Intra-class correlation coefficient values 

Raw data quantitative air trapping measures were 
assessed using MLD and LAA% having attenuation 
below -850 HU or between -850 and -910 HU. The 
match between the two observers yielded good to excel- 
lent inter-reader reliability for each (Additional file 3: 
Table SI). The lowest inter-class correlation coefficient 
between junior and senior radiologists was 0.91 for the 
fourth slice EXP_85o-9io- 



Table 1 Demographic and spirometric characteristics of 
study subjects 



Variables 


Mean 


Std 


Age (y) 


72,03 


8,08 


Height (cm) 


162,14 


8,93 


Weight (l<g) 


67,05 


11,54 


BMI (Kg/m^) 


25,53 


3,49 


FVC (!) 


3,28 


0,89 


FVC (% predicted) 


123,10 


20,86 


FEVl (!) 


2,37 


0,64 


FEVl (% predicted) 


89,31 


17,02 


FEVl /FVC 


72,69 


6,82 


MEF (!) 


2,84 


1,14 


MEF (% predicted) 


110,09 


40,98 


TLC (!) 


5,85 


1,17 


TLC (% predicted) 


110,95 


15,21 


RV(I) 


2,58 


0,61 


RV (% predicted) 


117,13 


26,21 


FRC (!) 


3,28 


0,73 


FRC (% predicted) 


111,87 


22,64 


CV(I) 


2,91 


0,80 


dN2 (% slope) 


3,02 


2,24 



Note. Data are means ± standard deviations for quantitative variables. 
Percentages are expressed as a ratio of measured to predicted values. 
FVC -forced vital capacity, FEVl = forced expiratory volume in 1 second, 
MEF- mid forced expiratory flow, TLC -tota\ lung capacity, RV- residual 
volume, F/?C afunctional residual capacity, Ci/=closing volume, c/A/2 = phase III 
nitrogen slope. 
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Figure 1 Box plots of air trapping represented by E/I^ld 4 
levels and mean of E/Imld- Significant differences between slices 
were observed (p = 0.0002). 



Spatial Heterogeneity 

As shown in Figure 1, CT air trapping assessed by E/I 
MLD ratio obtained at all four levels was statistically dif- 
ferent (p = 0.02) but all the four expiratory mean lung 
density values remained significantly correlated with 
dN2 as demonstrated by the Spearman correlation coef- 
ficient (Additional file 4: Table S2). 

Comparison of CT indexes for small airway obstruction 
assessment 

43 of 89 subjects (48,3%) had dN2 value above the cut- 
off (2.5% N2/1). 

All CT indexes (E/I, (E-I)/I and E-I LAA% or MLD 
correlated with small airway obstruction, as assessed by 
dN2 (Additional file 4: Table S2). Figure 2 shows the 
connection between dN2 and the CT criteria generated 
by hierarchical dendrogram analysis. The cluster located 
the furthest from dN2 was E-Imld- The three branches 



the most closely related to dN2 in the cluster dendro- 
gram (E/Imld, (E-I)/Imld> E-I-850-9io) were then included 
for further ROC curve analysis to determine the criteria 
that yielded the highest combined sensitivity and specifi- 
city and thus the highest diagnostic value. The resulting 
ROC curves are shown Figures 3, 4 and 5. E/Imld Pro- 
vided the best area under the curve (AUC) value (0.714 
confidence interval: 0.606, 0.822) using dN2 as a re- 
ference. Similar results were obtained for (E-I) /I. When 
analysis was performed on the second slice individually, 
AUC for E/I and (E-I)/I was slightly higher (0.733; 
0.630 - 0.836). 

Discussion 

This study reports a clear relationship between spirome- 
trically gated CT air trapping quantification and physio- 
logical measurements of the small airways using SENT. 
The strength of the CT quantification technique is its 
capacity to express small airway obstruction reliably 
using standard, commercially available, software. Indeed, 
software-assisted CT air trapping quantification was 
feasible in our study for all subjects who achieved ad- 
equate CT acquisitions with a high consistency between 
observers even for inspiratory images, matched visually 
by using anatomical landmarks. Expiratory to inspiratory 
ratio and the difference between these values were calcu- 
lated to characterize CT air trapping in order to over- 
come error related to decrease in lung density due to 
alveolar enlargement likely associated with age [26-28]. 
E/Imld and (E-I) /I mld exhibited the best performance 
characteristics to assess small airway obstruction. O. 
Mets and colleagues recently pointed out that E/Imld 
was related to PFTs [29] . They showed a relationship be- 
tween residual volume/Total lung volume capacity (RV/ 
TLC) and E/I mld in a population of current and former 
heavy smokers in a lung cancer screening cohort. None- 
theless, RV/TLC reflects static lung hyperinflation and 



Cluster Dendrogram 



Figure 2 Hierarchical dendrogram analysis to evaluate distance and similarity between dN2 and CT air trapping criteria. 
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0.6 0.4 
Specificity 

Figure 3 Receiver-operating characteristic (ROC) curves show 
the diagnostic performance of E/l using dN2 as a functional 
test of airway obstruction (dN2 cut-off = 2.5% N2/I). The best 
threshold is reported, alongside the sensitivity and specificity: 
E/l = 0.89 (0.70; 0.73) - E/l _85o-io24 = 0.32 (0.63; 0.67) - 
E/l -850-910 = 0.29 (0.79; 0.61). 




0.6 0.4 
Specificity 

Figure 5 Receiver-operating characteristic (ROC) curves show 
the diagnostic performance of (E-l)/l using dN2 as a functional 
test of airway obstruction (dN2 cut-off = 2.5% N2/I). The best 
threshold is reported, alongside the sensitivity and specificity: 
(E-l)/l MLD = -0.1 1 (0.70; 0.73) - (E-l)/l -850-1024 = -0.68 (0.63; 0.67) - 
(E-l)/l -850-910 = -0.70 (0.79; 0.61). 
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Figure 4 Receiver-operating characteristic (ROC) curves show 
the diagnostic performance of E-l using dN2 as a functional 
test of airway obstruction (dN2 cut-off = 2.5% N2/I). The best 
threshold is reported, alongside the sensitivity and specificity: 
E-l MLD = 90.01 (0.79; 0.63) - E-l -850-1024 = -34.24 (0.65; 0.65) - 
E-l -850-910 = -47.23 (0.81; 0.59). 



the value may be influenced either by the presence of 
small airway obstruction or by emphysema. Further- 
more, discrimination is possible only in severe cases of 
hyperinflation. Only 38 of the 427 heavy smokers sub- 
jects presented abnormal RV/TLC in the Mets study. 
There is currently no functional test allowing perfect 
non- invasive assessment of small airways. Nonetheless 
inert gas washout tests such as SENT provide infor- 
mation on ventilation inhomogeneity that allows sensi- 
tive detection of small airway obstruction even when 
RV/TLC or other standard spirometry values are normal 
[30-33]. Our sample population aged 60 and older and 
the selection criteria used offered a well-defined popu- 
lation of subjects with a range of small airway obstruc- 
tion that may be attributed to normal ageing [34,35]. 
Furthermore, people included in this study achieved 
reproducible maneuvers that are a key point for SENT 
validity [6]. 

The correlation found in our study between CT and 
SENT showed that CT air trapping can detect even mild 
to moderate small airway obstruction as seen in our sub- 
jects who were all asymptomatic but nevertheless having 
physiological modifications. The aim was to validate CT 
air trapping quantification from a homogeneous model 
that could be further applied to various chronic obstruc- 
tion situations such as COPD or asthma. 

Our study has several potential limitations. First, we 
did not used whole expiratory volume also described for 
air trapping assessment but only a single set of four 
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expiratory slices [19,36]. We did this because of the 
difficulty for subjects to sustain breath-hold at low lung 
volumes and furthermore in keeping with the recom- 
mendation made by the ethics committee to limit the 
radiation dose as well as the "As Low As Reasonably 
Achievable" principle. We cannot know from this study 
whether analysis on the whole lung volume would 
have improved the correlations, nonetheless our results 
showed similar significant air trapping heterogeneity be- 
tween the upper, middle, and lower regions to that found 
by Bankier and colleagues in their work on a lung trans- 
plant population [37]. Regardless of slice level, each 
remained significantly associated with dN2. Even if AUG 
was slightly higher for E/I mld in slice 2 (i.e. through the 
carina), we did not recommend the use of only one ex- 
piratory slice as its index could not be applied to path- 
ologies such as asthma, which have been proven to have 
widely heterogeneous spread [38]. 

Second, although use of a spirometer to control lung 
volume during acquisition is not standard practice, in 
this study it provided proof of inspiration and expiration 
volume. 

Another potential limitation is the lack of emphysema- 
tous lesions in our population. The elective site of meas- 
urement should be adapted in COPD patients with such 
lesions as stated by Matsuoka [39]. 

Conclusions 

Our study demonstrated that in asymptomatic subjects 
with small airway obstruction, software-assisted CT air 
trapping quantification correlated with SENT. E/I and 
(E-I)/I MLD are equally representative of small airway ob- 
struction and could be used as a small airway obstruc- 
tion biomarker. 

Additional files 
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